Noise-assisted spike propagation in myelinated neurons 
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We consider noise-assisted spike propagation in myelinated axons within a multi-compartment 
stochastic Hodgkin-Huxley model. The noise originates from a finite number of ion channels in each 
node of Ranvier. For the subthreshold internodal electric coupling, we show that (i) intrinsic noise 
removes the sharply defined threshold for spike propagation from node to node, and (ii) there exists 
an optimum number of ion channels which allows for the most efficient signal propagation and it 
corresponds to the actual physiological values. 
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I. INTRODUCTION 

In many vertebrates, the propagation of nerve im- 
pulses is mediated by myelinated axons. These nerve 
fibers are composed of active zones, the so-called nodes 
of Ranvier, where ion channels are accumulated, and 
passive fragments which are electrically isolated (myeli- 
nated) from the surrounding electrolyte solution. The 
myelin sheath presents a crucial evolutionary innovation. 
In a process called saltatory conduction, the neural 
impulse propagates from node to node more rapidly 
than it propagates in unmyelinated axons of equal 
diameter [l|, 0] . Myclination not only helped Nature to 
increase the signal propagation speed but also to dras- 
tically reduce the metabolic cost of neural computation. 
Indeed, with a single neural impulse (of order of 1 ms), 
each sodium channel transfers up to 10* sodium ions into 
the cell which then should be pumped back to restore 
the appropriate steady-state electrochemical potential. 
To transfer three sodium ions, the corresponding Na-K 
pump hydrolyzes one ATP molecule [2], which thus 
yields the estimate of about 3 • 10'^ ATP molecules per 
ion channel per neural spike. A too large number of 
channels imposes an extremely high metabolic load (the 
brain of the reader is just now consuming about 10% of 
the body's metabolic budget, which, per one kilogram of 
mass, is more than the muscles use when active). More 
elaborate estimates confirm that even a small cortical 
cell needs, in a long_£un, at least 10^ ATP molecules per 
one neural spike 



neurons tend to minimize the number of ion channels 
because of energetic costs. On the other hand, too 
few channels may trigger random, parasitic spikes or 
induce spike suppression by strong fiuctuations, making 
the signal transmission too noisy and thus unreliable. 
The questions we investigate in this paper within a 
simplified stochastic model of signal transmission in a 
myelinated neuron, are as follows: can channel noise 
soften the propagation threshold and facilitate the signal 
propagation which would not occur in the deterministic 
case? Does an optimum size of the channel cluster 
exist, which in turn yields a most efficient, noise- assisted 
propagation? 

The statement of this objective shares features similar 
in spirit to noise supported wave propagation in subex- 
citablc media (see Ref. @ and references therein) and 
noise in excitable spatiotcmporal systems (see Ref. 0). 
The infiuence of intrinsic noise on membrane dynamics 
was studied in the context of intrinsic stochastic reso- 
nance, see in Refs. and synchronization of ion chan- 
nel clusters [lo| . A new aspect emerging from the present 
study is a possible interpretation of the neural spike 
transmission as a delayed synchronization phenomenon 
occurring in the chain of active elements 



II. THE MODEL 



A. Compartment modelling for myelinated axons 



n pnenomei 



In myelinated neurons, yet another problem emerges: 
there is a threshold present for the electrical coupling 
between the nodes of Ranvier. Deterministic cable 
equation models predict that if the internodal distance is 
too large, the coupling becomes too small and the signal 
propagation consequently fails 0- However, due to a 
finite number of channels (of the order of 10*), intrinsic 
noise is inevitably present in the nodes of Ranvier Q. 
The deterministic models can only mimic the behavior 
of a very large number of ion channels possessing a 
negligible intrinsic noise intensity; in contrast, real 



As an archetypical model for signal transmission along 
a myelinated axon of a neuronal cell we consider a 
compartmcntal stochastic Hodgkin-Huxley model (see 
Sec. Ill BlIII C|) . In contrast to unmyelinated axons, where 
ion channels are homogeneously distributed within the 
membrane, the myelinated axon consists of alternating 
sections where the ion channels are densely accumulated 
(nodes of Ranvier), and regions with very low ion chan- 
nel density, encased in multiple layers of a highly resistive 
lipid sheath called myelin. The nodes of Ranvier play the 
role of active compartments in our model, while the elec- 
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FIG. 1: Sketch of a myelinated axon; The axonal cell membrane in the nodes of Ranvier (depicted in light-grey) contains a 
high ion channel concentration. The remaining segments are wrapped in the myelin sheath (dark-grey). Spike propagation 
along the axon occurs in a saltatory manner. 



trically neutral myelinated segments constitute passive 
compartments (see Fig. [T]). The typical intcrnodal dis- 
tance L is about 1-2 mm while the length A of the nodes 
of Ranvier is in the micrometer range [l| . 

The generalization of the original Hodgkin-Huxley 
model is related to the influence of the channel noise 
which results from the randomness of the ion channel 
gating. The membrane potential at each particular node 
of Ranvier is assumed to be constant across the whole 
node area and is characterized hy Vi, i = 0,1,2, ....N — 1 
where A'' is the total number of nodes [l|, Q . 

Supposing total electrical neutrality of the passive re- 
gions, the electrical properties of a myelinated axon are 
modeled by a linear chain of diffusively coupled active 
elements. The dynamics of the membrane potential is 
given by: 



axon as well as on the resistivity of the extracellular and 
intracellular medium, and serves as a control parame- 
ter in our studies It depends also on the ratio of 
the node's diameter d and its length A. Typically, the 
node's diameter is two orders of magnitude smaller than 
its length 0]. 

Note that the Eq. ([T]) just presents the Kirchhoff law 
for an electrical circuit made of membrane capacitors, 
variable nonlinear membrane conductances and intcrn- 
odal conductances k. 



B. Hodgkin-Huxley modeling 

According to the standard Hodgkin-Huxley model [if 
the ionic membrane current reads: 



C — Vi — /i.ionic(^, t) + Ii.cxt{t) 

at 

+ /,,inter, for I = 0, 1 , 2, . . . .TV - 1 , (1) 

where C denotes the axonal membrane capacity per unit 
area. Ii^ionic{Vi,t) is the ionic membrane current (per 
unit area) within the ith node of Ranvier. I i, ionic {Vi,t) 
depends on time t and on the membrane potential Vi in 
the specific node i only, and is described by a neuronal 
membrane model which we subsequently treat within a 
coupled stochastic Hodgkin-Huxley setup). /j;,oxt(i) de- 
scribes an external current per unit area applied on the 
ith node. In our setup wc apply the stimulus, i.e. a con- 
stant current, on the first node only, bringing it into a 
dynamical regime of periodic firing [14] . If the coupling 
between nodes is sufficiently strong, the action poten- 
tials may start propagating along the axon. The coupling 
to the next-neighboring nodes of Ranvier is achieved by 
inter-nodal currents (per unit area) /i, inter, given as: 

f K iV+i - Vi) for I = , 

4intcr = < -^,) for^ = iV-l, (2) 

Ik {Vi-i — 2Vi + Vi+i) elsewhere. 

The coupling parameter n depends, among other things, 
on the length L of the intcrnodal passive segment of the 



/^(t^^) =-GkK;) (V-Ek) 

- GNa(m„ h,) {V - ^Na) - Gl(F, - Ei^) , (3) 

with the potassium and sodium conductances per unit 
area given by: 

= n\, G^,{m„h,) = g^^ mlh, . (4) 

In Eq. (121), Vi denotes the membrane potential at the ith 
node of Ranvier. Furthermore, i?Na, E^ and i?L arc the 
reversal potentials for the potassium, sodium and leakage 
currents, correspondingly. The leakage conductance per 
unit area Gl is assumed to be constant. The parame- 
ters g^'^^ and g^^^ denote the maximum potassium and 
sodium conductances per unit area, when all ion channels 
within the corresponding node are open. The values of 
these parameters are collected in the Table U [l^ . Note 
that in the nodal membrane the conductances of open ion 
channels are supposed to be Ohmic-like, i.e. the nonlin- 
earity derives from their gating behavior only. Moreover, 
formulating the Eqs. ([3]) and we implicitly assumed 
for simplicity that the different axonal nodes are kineti- 
cally identical, i.e. the number of sodium and potassium 
ion channels is constant for each node. As a consequence, 
the maximum potassium and sodium conductances are 
identical constants for every node of Ranvier. 
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The gating variables rii, rui and hi, cf. Eqs. ([3]) and 
dH) , describe the probabihties of opening the gates of the 
specific ion channels in the it\\ node upon the action of 
activation and inactivation particles, h is the probability 
that the 1 inactivation particle has not caused the Na 
gate to close, m is the probability that 1 of the 3 required 
activation particles has contributed to the activation of 
the Na gate. Similarly, n is the probability of the K 
gate activation by 1 of the 4 required activation particles. 
Assuming gate independence, the factors and hi 
are the mean portions of the open ion channels within a 
membrane patch. The dynamics of the gating variables 
are determined by voltage-dependent opening and closing 
rates ax{V) and (3x{V) [x = m, h, n) taken at T = 6.3 °C. 
They depend on the local membrane potential V and read 
(with numbers given in units of [mT/]) (isl. [l7j: 



am{V) = 0.1 



y + 40 



ah{V) 
an{V) 

Pn{V) 



l-exp{-(y + 40)/10} ' 
4 exp{-(V^ + 65)/18} , 
0.07 exp{- (\/ + 65)/20} , 
1 

l + exp{-(V^ + 35)/10} ' 
V + bb 
1 -exp{-(F + 55)/10} ' 
0.125 cxp{-(F + 65)/80} . 



(5a) 

(5b) 
(5c) 

(5d) 

(5c) 
(5f) 



Such an approximation is valid for very large numbers 
of ion channels and whenever fluctuations around their 
mean values are negligible. 



Stochastic generalization of the 
Hodgltin-Huxley model 



Because the size of the nodal membrane is finite there 
necessarily occur fluctuations of the number of open ion 
channels. This is due to the fact that ion channel gat- 
ing between open and closed state is random. Accord- 
ingly, for finite size membrane patches like the considered 
nodes of Ranvier, there are fluctuations of the membrane 
conductance which give rise to spontaneous action poten- 
tials [J, H, H, [H, [TqI l20l| and the references therein. 

The number of open gates undergoes a birth-and- 
death-like process. The corresponding master equation 
can readily be written down. The use of a Kramers-Moyal 
expansion in that equation results in a corresponding 
Fokker-Planck equation which provides a diffusion ap- 
proximation to the discrete dynamics [l^, ll^l . The cor- 
responding multiplicative noise Langevin equations pll | 
then read: 



d 
dt 



Xi = ax{Vi) (1 - Xi) - PxiVi) Xi +£,i.xit) , (7a) 



The dynamics of the mean fractions of open gates reduces 
in the standard Hodgkin-Huxley model to relaxation dy- 
namics: 



—Xi = axiVt) (1 



Xi) - (3xiVi) Xi, x^m,h,n. (6) 



with X = m,h,n. Here, S,i.x{t) are independent Gaus- 
sian white noise sources of vanishing mean and vanish- 
ing cross-correlations. For an excitable node with the 
nodal membrane size A the non-vanishing noise correla- 
tions take the following form: 



{^r At) knit')) 



ApNa 
1 

ApNa 
1 

ApK 



[amm (1 - m,) + f3AV^) m,] S{t - t') , 
[aHm{l~hi) + hmhi\ 6{t~t'), 
[a„(yO (1 - Ui) + /?„(X^,) n,] 5{t - t') , 



(7b) 
(7c) 
(7d) 



with the ion channel densities being pNa and pK (see Ta- 
ble [l}. The stochastic Eq. dTaj) replaces Eq. dH) [111. 
Remarkably, the nodal membrane size A, which is the 
same for all nodes, solely influences the intrinsic noise 
strength. Note that the correlations of the stochastic 



forces in these Langevin equations contain the corre- 
sponding state-dependent variables and, being an ap- 
proximation to the full master equation dynamics, thus 
should be interpreted in the Ito sense [23| . 



III. OPTIMIZATION OF THE SIGNAL 
PROPAGATION 

We next study the model of a myelinated axon in which 
one node of Ranvier is continuously excited to sponta- 



neous spiking by an external supra-threshold current, 
whereas other nodes are initially in the resting state and 
the inter-nodal conductance k is too low for the spikes to 
propagate to the neighboring nodes in the deterministic 
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stochastic Hodgkin-Huxley model 
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Axonal chain model 
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Simulation time 


T 




3-10^ 


ms 


Simulation time step 


dt 




0.002 


ms 


Initial values for each node {i = 0, .., A'^) 


Voltage 


V^ 




-59.9 


mV 


Inactivation probability for the Na gate 


hi 




0.414 




Activation probability for the Na gate 


rrii 




0.095 




Activation probability for the K gate. 


rii 




0.398 





TABLE I: Model and simulation parameters 



model (notice the arrow in Fig. [2]). The fluctuations of 
the potential caused by the randomly opening and clos- 
ing channels can be large enough to help the spike to 
overcome the large inter-nodal resistance and to make 
the saltatory conduction between nodes possible. 

Deterministic limit 

We start out with the deterministic situation using a 
finite number of Ranvier nodes, i.e. the channel noise 
level is set to zero, which is formally achieved in the limit 
A oo. In order to provide some insightful explanation 
on the dependence of saltatory spike propagation on the 
coupling parameter, we numerically integrated the dy- 
namical system, see. Eqs. (UJ-® with the parameters 
given in Tab. |T1 In order to achieve equilibration along 
the whole chain we do not apply the constant current on 
the first node /o,ext from the beginning. Instead, we use 
the following protocol: (i) we initially allow for equilibra- 
tion of every individual node by integrating over 100 ms 
with K = 0; (ii) in a second step we switch on the cou- 
pling; (iii) finally after integrating over 150 ms, i.e. under 
stationary conditions, we apply a constant current stim- 
ulus of /o,cxt = 12/xA/cm^ on the first node only. Due 
to the supra-threshold stimulus [H] on the first node, the 
dynamics of the membrane potential of the initial node 
Vq exhibits a limit cycle, resulting in the periodic gen- 
eration of action potentials. The numerically obtained 
spike trains Vi for i = 0, .., 9 allow for studying the prop- 
agation of action potentials along the linearly coupled 
chain. Using a threshold value of Vth = 20 mV enables 



the detection of an action potential (also called spike or 
firing event) in the particular spike train of the individual 
nodes. Since we are interested in the successful transmis- 
sion along the whole chain, the spiking at the terminal 
node is analyzed and the number of spikes at the termi- 
nal node is related to the number of initiated spikes at 
the initial node. This defines the transmission reliability. 

For the chosen coupling strength k, the excitation 
propagates along the chain of given = 10 nodes of 
Ranvier. This particular choice of the number of nodes 
is somewhat arbitrary, although physiologically realistic. 
Note also, that these obtained results are robust and are 
typical for larger node numbers as well (not shown) . Ob- 
viously, for too small coupling parameters no excitation 
proceeds to the neighboring node "1" and, at the end, 
will not reach the terminal node "9" . In this case, we ob- 
serve that spike propagation fails. In the opposite limit, 
i.e. for large coupling constants /t, each action potential 
travels along the chain and arrives at the terminal node 
"9" . Then the spike transmission efficiency of 100% is 
achieved. In this situation, the dynamics of the terminal 
node is synchronized with a delay with the initial node, 
i.e. Vg{t) = Vo{t — t), where the delay time r accounts 
for the finite propagation speed. Most interestingly, how- 
ever, there occurs no sharp transition between those two 
regimes of 0% and 100% signal transmission (see Fig.[2|). 
For ftci = 0.0665 mS/cm2 and = 0.1360 mS/cm^, 
there is an intermediate regime KcI < k < Kc2 where 
discrete transmission patterns k : I occur. Here, k is the 
number of spikes generated at the initial node, while I 
correspond to the number of spikes transmitted to the 
final node. In principle, the coupling parameter k could 
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FIG. 3: The simulated firing events vs. time, cf. Eqs. JTJ - 
(H}, (|7ap -(f7l). of the noisy membrane potentials at three differ- 
ent nodes are depicted for k = 0.065 mS/cm^ and the nodal 
area A = 10* /xm^. The constant current stimulus /o,cx keeps 
the initial node "0" periodically firing. Noise-assisted spike 
propagation occurs only for some spike and some spikes do 
not reach the terminal node "9" . 
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FIG. 2: Deterministic spike propagation in a myelinated axon 
containing A'" = 10 Ranvier nodes (no channel noise) ; (a) The 
ratio of spikes generated in the node during an interval of 
1000 ms which propagate up to the node 9 (solid line), as a 
function of the value of the inter-nodal coupling strength (i.e. 
the conductance) k. The arrow indicates the subthreshold 
value of K = 0.065 mS/cm^ chosen for further simulations 
of the influence of channel noise on spike propagation. The 
stepwise shape of the graph depicts the different transmission 
patterns: (b) 2:1 ; (c) 3:2 ; (d) 4:3 and (e) 1:1 . For example, 
for the ratio of 2:1, every second spike which is generated at 
n = propagates. (Small irregularities in step levels occur as 
a result of a finite counting statistics of propagating spikes.) 



be tuned in such a way that any rational transmission 
pattern k : I Yg, attained. This manifests as generahzed, 
delayed synchronization. 



Note that a similar effect shows up when one drives 
the Hodgkin Huxley system (i.e. a single node) with an 
ac-sinusoidal signal, where the ratio of spiking events to 
driving periods exhibits a rational number [24| . 



Influence of channel noise 

Next, we investigate how the nodal membrane size A, 
determining the channel noise intensity, influences the 
propagation of spikes along the axon. Towards this objec- 
tive we numerically integrated the Hnearly coupled chain 
model with the nodes treated by the stochastic Hodgkin- 
Huxley model, cf. Eqs. ^ - (P and Eqs. ((7a)) -((7)). 

Following the same protocol as for the deterministic 
case, we computed the probability for a generated spike 
in the initial node to be transduced to the terminal node. 
In our example the axon consists again of 10 nodes. For 
the coupling parameter, i.e. the inter-nodal conductance, 
we chose a sub-threshold value of k = 0.065 mS/cm^, 
i.e. below the critical value KcI for a deterministic spike 
propagation along the axon, see Fig. [21 It turns out that 
for coupling parameters slightly smaller than the critical 
value Kci the influence of the channel noise is most strik- 
ing. Due to the channel noise, we observe spike propaga- 
tion even for the sub-threshold coupling, i.e. for k < Kci- 
The intrinsic noise weakens the strict threshold and even 
for sub-threshold values of the coupling parameter, there 
is a no n- vanishing probability for spike propagation. In 
Fig. [3] the spike train for three different nodes is shown 
for K = 0.065 mS/cm^ and A = 10* ^m^, i.e. for an 
intermediate intrinsic noise level. The channel noise fa- 
cilitates the propagation of action potentials along the 
axonal chain. 

Fig. m depicts cases of spike propagation in the pres- 
ence of noise of different intensities for the chosen sub- 
threshold coupling parameter k. Apart from the noise fa- 
cilitated spike propagation one can identify another par- 
ticular phenomenon: The failure of spike propagation due 
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FIG. 4: Noisy spike propagation for a myelinated neuron with 
10 Ranvier nodes and for k — 0.065 mS/cm^. Single real- 
izations of noise-assisted spike propagation for varying nodal 
membrane size A. (a) A = 250 /xm^: Strong channel noise in- 
hibits strongly spike propagation so that propagation occurs 
only over a short distance, (b) A = 10* nm^: An intermedi- 
ate noise allows for the propagation over the longest distance, 
(c) ^ = 5 X lO** nm^: Weak channel noise only rarely allows 
to overcome the inter-nodal resistance. 



numbers of ion channels) the spike initiated at the first 
node more likely propagates to the next nodes. However, 
it is also likely that the spike collides with a particularly 
large fluctuation and thus becomes deleted. Tracking the 
behavior of a given node, one observes the skipping of 
some firing events [1^ , see Fig. [31 The propagation dis- 
tance is then quite short and the excitations rarely arrive 
at the terminal node, see Fig.[Ha). Note also, that for a 
high level of channel noise spontaneous spikes can occur 
at any node. These parasitic spikes are not triggered by 
a preceding spike in the neighboring nodes. Therefore, 
they deteriorate the information transfer along the axon 
and are not considered for the spiking statistics at the 
terminal node. 

Too weak noise (i.e. a large number of channels) does 
not allow for effective spike transmission to neighboring 
nodes. Only the very first excitation is transmitted over 
a larger number of nodes, while the propagation of the 
successive action potentials is rarely facilitated by noise, 
see Fig.lU^c). This specific enabling of propagation of the 
first spike only is due to the initial setup of the problem 
(initial conditions), because of the existence of a refrac- 
tory time period. 

For an intermediate channel noise intensity, however, 
the propagation distance is the longest, although spikes 
propagate less frequently than in the presence of an in- 
tensive noise, see Fig.UJ^b). Overall, one may expect that 
there is an optimum dose of internal noise for which the 
spike propagation along an axon is most probable. To 
validate this assertion, we determined the probability for 
an excitation stimulated in the initial node, to arrive at 
the terminal node being which node "9" in our case. The 
fraction of spikes arriving at a specific node depends on 
both, the distance which has to be covered and the noise 
level. With increasing distance, i.e. the node number, 
this fraction declines monotonically, cf. Fig. [U^a). The 
decay depends on the channel noise level, i.e. the nodal 
membrane size A. There is an intermediate noise level 
for which the decline is slowest. As a result, there exists 
an optimum dose of internal noise (in our example for the 
nodal membrane size near A w 3800/^m2) at which the 
signal transmission to the terminal node is most effective; 
i.e., the ratio of transmitted spikes assumes a maximum, 
cf. Fig. [5Jb). We emphasize that the occurrence of the 
optimal dose of intrinsic noise is robust under a change 
of the total number of nodes (not shown). Note that 
the observed maximum corresponds to a nodal number 
of sodium ion channels of iVNa = 2.28 • 10^ which corrob- 
orates with the physiological range of the number found 
for sodium ion channels in a node of Ranvier (experimen- 
tal studies report TVNa ~ 10^) 



IV. CONCLUSIONS 



The optimization of the signal propagation along a 
to a large channel noise level. For rather strong intrinsic myelinated axon of a finite size occurs as the result of 
noise (for small nodal membrane sizes A., yielding small the competition between the constructive and destructive 
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FIG. 5: Noise-assisted spike propagation for spikes initiated in tlie initial node ("0") out of 10 nodes of a linearly coupled 
chain: (a) The numerically obtained percentage of spike arrivals at the particular nodes and (b) its dependence on the nodal 
membrane size A. The chosen coupling strength k — 0.065 mS/cm^ is sub-threshold. There occurs an optimum noise intensity 
(corresponding to a nodal membrane size A ~ 3800/im^) at which the signal transmission along the whole chain is most efflcient. 



role of channel noise occurring in the nodes of Ranvier. 
On the one hand, when the number of channels in the 
node is very small, the strong fluctuations of the activa- 
tion potential facilitate a spike propagation among nodes 
which otherwise does not occur in the absence of chan- 
nel noise. On the other hand, at a high level of channel 
noise spikes cannot travel over a long distance because 
it is likely that a spike soon collides with another noise- 
induced spike leading to their mutual annihilation, or the 
noise can also suppress spike generation. If the number 
of channels in nodes is large, the fluctuations are too 
weak for a spike to overcome the internodal resistance 
and subsequently to propagate. For a certain interme- 
diate number of channels in the node of Ranvier node, 
however, the intrinsic noise is sufficiently strong to in- 
duce the saltatory conduction while being still sufficiently 
weak not to deteriorate the spiking behavior. An optimal 
nodal membrane size can be identified, for which the sig- 



nal becomes most efficiently transmitted over a certain 
axonal distance. Moreover, the corresponding number of 
sodium ion channels corresponds to actual physiological 
values. This feature is quit e in spirit of the stochastic 
resonance phenomenon |28l . [29j with an intrinsic noise 
source d, Q. One may therefore speculate, whether Na- 
ture adopted this optimization method to balance the 
signal transmission efficiency and the metabolic cost. 
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